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Abstract 

The motion of I = antibound poles of the S-matrix with varying potential strength is calculated 
in a cutoff Woods-Saxon (WS) potential and in the Salamon- Vertse (SV) potential, which goes 
to zero smoothly at a finite distance. The pole position of the antibound states as well as of the 
resonances depend on the cutoff radius, especially for higher node numbers. The starting points (at 
potential zero) of the pole trajectories correlate well with the range of the potential. The normalized 
antibound radial wave functions on the imaginary fc-axis below and above the coalescence point 
have been found to be real and imaginary, respectively. 

PACS numbers: 21.10.Pc,21.30.-x,21.60.Cs 
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I. INTRODUCTION 



Nuclear states are most often described in terms of single-particle (s.p.) bases generated 
by a spherical potential, mostly of Woods-Saxon (WS) type. Bound and discrete unbound 
s.p. states all obey the outgoing-wave boundary condition, which is u(r, k) ~ exp(ikr) when 
both the charge and the angular momentum / are 0. The general solution behaves like 
exp(— ikr) — Sexp(ikr), where S, a function of the energy E or the wave number k, is called 
the S- "matrix" . Where the outgoing boundary condition is satisfied, the S-matrix has poles. 
The bound-state poles belong to E < or imaginary wave number with Imk = 7 > 0. The 
resonance poles belong to complex E and k, with k = ±k — vy (k, 7 > 0). For antibound 
(virtual) states, E < 0, k — — ij (7 > 0). 

A WS basis is only complete if, in addition to bound states, it contains continuum 
scattering states and/or resonances and/or antibound states [l|, Q|. The completeness is 
understood with respect to a generalized scalar product. The resonance states, which have 
definite intuitive meanings, have proved to be very useful in describing weakly bound or 
unbound states of nuclei [2j, unlike antibound states, whose exponential tail, exp(7r), looks 
unphysical. However, the inclusion of an antibound state of 10 Li 3| in the description of 
n ' 12 Li was found meaningful ^4 - 7 ] . This shows that antibound states and the corresponding 
S-matrix poles ( "antibound poles" ) do deserve some attention. 

nn 

As an extension of recent studies [8|, [9j of the dependence of the S-matrix poles on the 
tail behavior of the potential, we now study antibound poles. The nuclear potential should 
in principle have an exponentially decreasing tail, like the folding of the nuclear matter 
density with the one-pion exchange force. The standard WS potential obeys this criterion, 
but it can only be treated properly in anal ytic al calculations, and analytical solution to the 
Schrodinger equation with a WS potential [lCj only exists for angular momentum zero. The 
matter is that in numerically solving the problem with a prescribed boundary condition 
the solution has to be matched, at a finite distance, to the solution with potential zero 
(asymptotic solution), and the matching amounts to cutting off the tail of the potential. 
The error committed in this way is usually believed to be small, but in a recent paper it 
was shown that, for broad resonances, the poles in a cutoff WS potential strongly depend 
on the value of the cutoff radius 0, 9]. 

In this work we examine the effect of the cutoff on the WS potential, and compare its 



2 



behavior with a potential that goes to zero smoothly and is exactly zero beyond a point. 
This "Salamon-Vertse (SV) potential" contains as many parameters as a cutoff WS potential 
and its shape is similar except for its tail. The tail of the SV potential can only conform 
to that of the WS at the expense of the inner region. Conformity in a longer section can 
be achieved with more parameters. This paper is only concerned with pointing out where 
problems might appear because of the cutoff. 

Unlike in former studies of antibound states we are aware of , we now explore the 

wave functions as well. We limit our attention to I = since antibound states may only 
play some role for s-states. 



II. POTENTIALS 



In solving the radial Schrodinger equation, a numerically calculated inner solution has 
to be matched at a distance r = r a to the solution of the asymptotic equation, and that 
yields the S'-matrix. This procedure is tantamount to cutting off the WS potential at 
r = -R max < r a . The potentials will be given in a form that expresses that they are exactly 
zero beyond a point, i.e., they are of finite range in a strict sense. The cutoff WS potential 
is thus 

V ws (r) = K,/ WS (r), (1) 

with 

/ ws (r)= |-( 1 + e ^) 1 > if ^<^ (2) 
[ , if r > R max . 

In the resonance region the pole trajectories obtained by varying Vq do depend on the cutoff 
radius i? max 

The SV potential becomes zero beyond a finite value r > po such that all its derivatives 
are also zero. Thus the potential is differentiable in the whole domain r £ [0, oo), in contrast 
with the cutoff WS potential, which has a discontinuity at the cut. 

To follow Eq. (CQ), we write the SV potential as 

V SV (r) = V f sv (r), (3) 

where 

f V (r) = f sv (r, Cl ,po lPl ) = f P0 (r) - Cl f'(r), (4) 
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with 

fp{r) = < ' <- r >) 




fir) = J ^ (6) 



The range parameters po an d Pi are chosen as po > Pi, thus the potential in Eq. ([3]) vanishes 
at Pq. To make the SV potential conform to the WS potential, we fit its three parameters, 
p , pi and d (ci > 0), to the WS form / ws (r) [9j. 

To have several antibound poles in the same potential, we choose the neutron potential 
to represent a heavy nucleus, 208 Pb. The values R = 1.27 x 208 1 / 3 fm = 7.525 fm, a = 0.7 



fm were adopted [12|, with -R max = 15 fm. The SV parameters giving the best fit to the WS 



shape are: Ci = 0.997, p = 10.963 fm and pi = 8.328 fm 



III. WAVE FUNCTIONS 

Let us sketch briefly how the pole solutions of the radial equation are calculated. For 
I = the radial equation is 

d -^- + [k 2 - U(r)]u(r, k) = , (7) 

where U(r) = (2p/h 2 )V(r). We introduce an intermediate distance R{ m , where the internal 
("left") and external ("right") solutions are to be matched. The left solution is defined in 
the interval r G [0, i?; m ] such that 

ob left (r, k) 



Wleft(0,fc) = 0, 



dr 



= u'(0,k) = l. (8) 



r=0 

The right solution is defined in the interval r e [-Rim, ?"a]j where r a is in the asymptotic region 
(r a > -R m ax and r a > p ), so that the solution satisfy the boundary condition 



u 



right 



r a ,fc) = e^ {u' Tight (r a ,k) = ike ikr «). (9) 



We integrate Eq. (J7|) numerically starting from the origin up to _Ri m and from r a down to 
.Rim- The eigenvalue or pole position is defined as the k value for which the left and right 
logarithmic derivatives 

j (h \ _ u Mt( R im,k) u' Tight (R iia , k) 

Mleft { -Kim ,K) Wright [ -Kim , K ) 



are equal: 

Lieft(kj) - L light (kj) = 0, (11) 

where, for antibound states, kj = —i^j (jj > 0), with j denoting the sequence number of 
the state. If we introduce the matching factor ai e ft = w r ight(-Rim, kj)/u\ e { t (Ri m , kj) of the left 
solution, the eigensolution 

. , , I aieftMieft(r, kj) , if r < R im 
v{r, kj) = < (12) 

[ Wright (r, kj) , if r > R ini 

obtained in this way is well matched but not normalized. 

Since f r °° exp(2ikr)dr = oo for k = — 27 with 7 > 0, the norm of an antibound state is 
infinity in the normal sense. By truncating the norm integral / °° v 2 (r, kj)dr at r = r a , the 
result will depend on r a through a term (2ikj)~ l exp(2ikjr a ), which has to be eliminated. For 
resonance states this term can be eliminated either by using the prescription of Hokkyo [13I , 
or by rotating the integration path of f°° exp(2ikr)dr onto the complex r-plane to the extent 
that the primitive function go to zero in infinity [14], which results in — (2'ik)^ 1 exp(2ikr a ) 
for the integral, and cancels the spurious dependence on r a resulting from J Q ra t> 2 (r, kj)dr. 
This rotation of the integration path provides a sound generalization for the scalar product 
involving Gamow resonances [l|, and makes it possible to construct complete sets involving 
resonance states. The same prescription also sets the tail term of the norm integral of 
an antibound state to — (2ik)~ 1 exp(2zfcr a ) if a more radical rotation (by an angle > it) is 
applied, and the results with this formula are meaningful [ijj]. It is this prescription that 
allows the inclusion of antibound states in complete sets of states 161 ] . 

With this, the square of the norm of v(r, kj) is 

N 2 = r v 2 {r,kj)dr-C(r a ,kj), (13) 
Jo 



where 



2ikr 

C(r,k) = —. (14) 



The antibound wave function normalized to 1 is thus 

u(r,kj) = —v(r,kj). (15) 

For kj = —i^ij (7j > 0), the term C(r a ,kj) is positive, just as the first term in Eq. (fi3|) . 
Thus iV 2 may be either positive or negative, a fortiori N as well as u(r, kj) may be real or 



imaginary. Since the radial wave function u enters in the norm integral as u 2 , and u must 
not be complex conjugated in any matrix elements the imaginary wave function causes 
a strange behavior [4]. 

The pole positions kj and the corresponding normalized radial wave functions were cal- 
culated by a modified version of the computer code GAMOW [17| . The accuracy of the 
calculation was checked by a more accurate program ANTI [18] using Ixaru's CP method 
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IV. NUMERICAL RESULTS 

A. Qualitative behavior of antibound poles 

Figures Q] and [2] show the imaginary part of the pole wave number A; as a function of 
the potential depth for the WS and for the SV potential, respectively. For bound and 
antibound states Re (k) = 0. For a very shallow potential, there is just one antibound state, 
with node number n = 0. With the attraction increased, the pole passes through the origin 
at Vo = Vo,o, and the system becomes bound. The Im (k) versus Vb curves belonging to the 
other poles look like parabolas with horizontal axes. The bound states become antibound 
as the potential depth is decreased to Vb jn , and meet another antibound pole at V n . What 
happens beyond their coalescence can only be depicted on the complex fc-plane (Fig. [2D- 



The two poles part the Im (fc)-axis perpendicularly in opposite directions 20] . 

We thus see that, while the bound state poles all move upwards along the imaginary fc-axis 
when the potential is deepened, some antibound states behave conversely. The energy shift 
caused by a perturbation SV f(r) can be estimated by 5E = u 2 (r,kj)5V f(r)dr. The 
sign of SE with respect to that of 5Vof(r) depends on whether u(r, kj) is real or imaginary. 
By looking at the Vo dependence of the pole, we can unambiguously infer that the wave 
function is imaginary on the upper branches of the parabolas, and they are real below. The 
single n = antibound-state wave function is imaginary. 

In Fig. H] we show the radial wave functions of some normalized antibound states in WS 
potentials. The antibound states that belong to the same node number in two different 
branches of the parabola seem to be non-orthogonal to each other although they are gen- 
erated, pairwise, by the same potential. That is, however, just an appearance. In fact, the 
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FIG. 1: Imaginary part of the pole wave number as a function of the depth of the WS potential. 
For bound and unbound states Im (k) > and Im (k) < 0, respectively. 
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FIG. 2: Imaginary part of the pole wave number as a function of the depth of the SV potential. 

tail region of the overlap integral cancels the the contribution of the inner region. The j = 2 
and j = 3 antibound states are orthogonal to each other, and so are the?' = 4 and j = 5 
states. Thus, pairwise, they may be included in complete sets of states |l] simultaneously. 
(The antibound states of different node numbers are, of course, orthogonal to each other if, 
unlike in Fig. HJ they are produced by the same potential.) 

If we have a centrifugal or Coulomb barrier, the picture is different in that the bound-state 
poles meet the antibound poles at the origin, and bifurcate there into a pair of resonance 
poles. In the (Vo,Im (k)) plane this corresponds to parabolas whose apices are at the origin. 
When the potential bottom is lifted, the antibound poles approach the origin monotonously 
from below, thus their normalized wave function is real throughout. 
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FIG. 3: Trajectories of the two n — 1,1 — poles in the WS potential with Vq varied (-R m ax — 

15 

fm). 

B. Quantitative observations 

A numerically most sensitive quantity is the apex V n of the parabolas in Fig. [IJ and 
that was used for testing the _R max - dependence for the WS potential (Tabled]). We see that 
for n = 1,2,3 the V n values are practically independent of -R max - The largest variation 
is in V3, most probably due to the enhancement of the error in the numerical solution of 



the differential equation as discussed in Ref. 21] . The k- values of the apices are somewhat 



more sensitive to R max , and the sensitivity gets more pronounced for higher n. (We will 




r [fm] 



FIG. 4: Normalized radial wave functions of antibound states in WS potentials. The n = 1 
(j = 2, 3) and n = 2 (j = 4, 5) states were produced by Vo = 6.9692 and 19.5 MeV, respectively. 
The wave numbers kj (in fm" 1 ) are k 2 = -0.183 i, k 3 = -0.188 i, k A = -0.410 i, k 5 = -0.422 i. 
The functions with j = 2,4 are imaginary, while those with j = 3, 5 are real. 



8 



return to this problem in discussing i? max - dependence of the pole trajectories, see Fig. [7] 
later.) The sensitivity to the potential shape has also been tested by comparing the values 
obtained for the potential strength Vo, n , which puts the pole at the threshold (Table ITT]). For 
WS, -R max = 15 fm was used, but it was ascertained that Vo jn is practically independent of 
-Rmax £ [15, 25] fm. The strengths for the two potential forms are very similar, which follows 
from the shapes being very similar. 

TABLE I: Well depths V n at the coalescence of the two antibound poles 



iW (fm) V l (MeV) V 2 (MeV) V 3 (MeV) 



WS 


15 


6.969 


18.995 


36.286 




16 


6.969 


18.992 


36.263 




18 


6.968 


18.989 


36.239 




20 


6.968 


18.989 


36.239 




25 


6.968 


18.989 


36.230 


sv 




6.978 


19.378 


37.347 



TABLE II: The values Vo jn setting the pole at the threshold 
Potential V 0fi (MeV) V ,i (MeV) F ,2 (MeV) V 0:3 (MeV) 
WS 0.897 7.727 20.562 39.122 

SV 0.893 7.634 20.519 39.072 



As we showed in Fig. [31 beyond the coalescence, the pair of antibound poles is transformed 
into a pair of decaying and capturing resonance poles. We show the trajectories of some of 
the / = decaying resonances in the complex fc-plane in Fig. [5j (The poles of the capturing 
resonances are the mirror images of the decaying ones with respect to the Im (fc)-axis.) 

The starting point of a trajectory is defined by the limit kj = liniy _,.o kj(Vo). For a 



potential of range R, an estimate for this limit is given by 



20] 



HIT 

Re(k n ) = — + 0(1). (16) 
For large n we can perhaps neglect the term 0(1). We approximate Vq = by 5 keV. 
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FIG. 5: Trajectories of the odd-ra, I = poles in the SV potential with Vo varied. The even-n 
trajectories are similar, except for n = 0, which runs along the imaginary fc-axis. 

In Fig. |5]one can see the trajectories of the I = 0, n = 1, 3, 5, 7 poles in the SV potential. 
Only the (anti)bound states have definite node numbers n, but resonances can also be 
characterized by the node number of the (anti)bound state that they correspond to. The 
real parts of the starting points are seen to be almost equidistant. Therefore, these Re (k n ) 
values can be fitted well by the straight line: Re (k n ) = a® + ain, with a slope a% = 0.32 
fm -1 , which implies R = 9.778 fm. 
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FIG. 6: Re (k n ) values of starting points of resonance trajectories for the WS potential cut off at 
-Rmax = 15 fm fitted with a straight line 

As for the WS potential, in Fig. E]one can see the starting points Re (k n ) as a function 
of n, and a straight line fitted to it. Although, for n < 4, the Re (k n ) values are somewhat 
erratic, the slope of the line, a\ = 0.212 fm -1 , provides R = 14.82 fm, in good agreement 
with -R max = 15 fm. To see the dependence of the R value deduced in this manner on 
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■■■■ n = 3 
-- n = 5 





-Rmax, we repeated the calculations for a set of -R max values chosen from the range typically 
used in practical calculations. The slope of the line was determined from five points with 
n = 4, . . . , 8. The results are given in Table UTTl The smaller -R max , the better is the agreement 
with R, and the better is Eq. ( 1T5|) satisfied. For larger -R max the rund-off errors of the 
numerical solution of the radial equation get larger. This fact forbids one to go substantially 
beyond i? max = 20 fm. 

TABLE III: Ranges R obtained from Eq. (|16p for different values of the cutoff radii i? max . The 
trajectory starting points Re (k n ) were fitted by a linear function of n, and R was calculated from 
its slope. The a values show the quality of the fit of the data to the straigt line. 



^ max (fm) R (fm) cr 



11 


10.93 


1.6 x 10~ 6 


14 


13.85 


2.0 x 10~ 5 


17 


16.89 


2.0 x 10~ 5 


20 


20.45 


6.4 x 10~ 4 



We examined the sensitivity of the pole trajectories to the cutoff radius, and in Fig. [7] 
we illustrate the results with the case of n = 7, which, for -R max = 15 fm, fits well into 
the straight line in Fig. |6j In view of the approximate independence of V n on -R max (see 
Table HJ, the results look surprising. We see that the n = 7 trajectory and, indeed, its point 
of intersection with the imaginary fc-axis, depend appreciably on the cutoff. While -R max is 
changed between 11 fm and 20 fm, the intersection of the trajectories with the Im (fc)-axis 
moves from — i0.40 fm -1 to — i0.61 fm -1 , with the potential depth to be set to 166 and 159 
MeV, respectively. Thus, similarly to the n = 1, 2, 3 cases, Vj is less sensitive to -R max than 
the pole positions. 

The stability of V n as a function of i? max can be understood, again, in a perturbative 
picture. The shift of a pole energy caused by changing the cutoff radius -R max from Ri to 
i?2 can be estimated to be AE = V ws (r)u 2 (kj,r)dr. Now, the tail of V ws (r) is small, 
but, for a resonance, u 2 (kj,r) ~ e 2lkr (for r £ [i?i,i?2]) is complex and may take large 
absolute values, and, correspondingly, the resonance poles may be shifted appreciably in 
the complex i?-plane as well as in the /c-plane. For antibound poles, however, the function 
u 2 (kj,r) is real, so that the pole can only be shifted along the imaginary fc-axis. In the 
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perturbative approximation the coalescence point of a resonance trajectory is thus shifted 
into the coalescence point of the shifted trajectory with unchanged V n , which suggests that 
V n need not be changed much when -R max is varied even in an accurate calculation. 

Looking at the trajectories in Fig. we see that the larger the value of -R max , the farther 
from the origin do they intersect with the Im (fc)-axis. Moreover, near vanishing potential, 
all trajectories start with a vertical section at a certain Re (kj) = K7. The larger the value 
of -R m ax, the smaller is Kj, and the inverse proportionality expressed by Eq. (JIB]) is borne 
out. 

It is interesting to compare this behavior with the case of the square-well potential ex- 
plored in Ref. [22j . For such a potential with radius R, the value of (3 n (R) = 7„i? (with 
— i% denoting the coalescence point) was found to be equal to 1, at least for low n values, 



independently both of R and of the node number n [TJJ. For a cutoff WS potential the 
corresponding (3 n = -R max 7n does depend on R max , and, for the n = 7 case shown in Fig. [7J 
can be approximated by a first order polynomial: (3 7 (R max ) = —5.05 + 0.864_R max . 
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FIG. 7: n = 7 resonance trajectories for the WS potential cut off at different i? max values. 



V. SUMMARY 

We can summarize the results as follows. 

The strange behavior of the antibound basis state found in Ref. |4J is explained by its 
normalized radial wave function being imaginary. Except for n — 0, the poles occur pairwise, 
and there is a range of potential depths in which there are two antibound states of the same 
node number: one below, and the other above the coalescence point. It has been shown 
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that the antibound states lying below the coalescence points are real, while those above are 
imaginary. This seems to be a general property of antibound wave functions. The antibound 
states may be included in an orthonormal basis. Numerical examples show that even those 
that belong to the same node number are orthogonal to each other. 

The pole belonging to node number n = is an exception; it starts (with an infinitesimally 
small attractive potential) as an antibound state, and becomes bound when the potential is 
deepened, without ever passing into the resonance region. The behavior of all other poles 
show similarity to the I > case the real parts of the starting points of the resonance 
trajectories (near potential zero) are inversely proportional to the potential range. For the 
WS potential, this range is to be identified with the cutoff radius. For the WS potential the 
pole trajectories, including the positions of the antibound states, depend on the cutoff radius, 
and the higher the node number, the stronger the dependence is. Thus, without discrediting 
the use of the WS potential in representing the nucleus in bound-state or scattering problems, 
this paper cautions against its indiscriminate use to represent broad resonances or antibound 
states. 
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